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Abstract. 

We consider the spin-1/2 Heisenberg chain with alternating spin 
exchange in the presence of additional modulation of exchange on 
odd bonds with period three. We study the ground state magnetic 
phase diagram of this hexamer spin chain in the limit of very strong 
antiferromagnetic (AF) exchange on odd bonds using the numerical 
Lanczos method and bosonization approach. In the limit of strong 
magnetic field commensurate with the dominating AF exchange, the 
model is mapped onto an effective XX Z Heisenberg chain in the 
presence of uniform and spatially modulated fields, which is studied 
using the standard continuum-limit bosonization approach. In absence of 
additional hexamer modulation, the model undergoes a quantum phase 
transition from a gapped phase into the only one gapless Liittinger liquid 
(LL) phase by increasing the magnetic field. In the presence of hexamer 
modulation, two new gapped phases are identified in the ground state 
at magnetization equal to ^ and | of the saturation value. These phases 
reveal themselves also in magnetization curve as plateaus at corresponding 
values of magnetization. As the result, the magnetic phase diagram of the 
hexamer chain shows seven different quantum phases, four gapped and 
three gapless and the system is characterized by six critical fields which 
mark quantum phase transitions between the ordered gapped and the LL 
gapless phases. 



PACS numbers: 75.10.Jm, 75.10.Pq 
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1. Introduction 

Plateaus observed in tlie zero-temperature magnetization curve of spin systems, usually 
display the quantum nature of this phenomenon. Formation of a plateau is connected 
with generation of a gap in the excitation spectrum, which can be in some senses 
regarded as a realization of generation the Haldane conjecture [1]. The quantum nature 
of a plateau formation mechanism was clearly shown in the seminal paper by Oshikawa, 
Yamanaka and Affleck [2], who proposed the condition for realization of a plateau at 
magnetization m as n{S — m) = integer, where S is the magnitude of the local spin, n is 
the number of spins in a translational unit cell of the ground state and m is normalized 
to 1 at saturation. 

A particular realization of such scenario appears in the one-dimensional (ID) space- 
modulated (alternating) quantum spin systems. The bond alternating Heisenberg spin- 
1/2 chains which are obtained by a space modulation in the exchange couplings represent 
one particular subclass of low- dimensional quantum magnets which pose interesting 
theoretical [3-17] and experimental [18-26] problems. The bond alternating spin-1/2 
chains have a gap in the spin excitation spectrum and reveal extremely rich quantum 
behaviors in the presence of an external magnetic field. By turning the magnetic field, 
the excitation gap reduces and reaches to zero at the first critical field. Simultaneously, 
the magnetization remains zero up to the first critical field. By more increasing of the 
magnetic field, system remains in gapless phase and at the second critical field, the 
gap re-opens and the saturation plateau is appeared in the magnetization curve. These 
models have only two plateaus at zero and saturation values of the magnetization. 

The subject of the space-modulated spin chains has grown in recent years and 
has been found that a space modulation in the exchange couplings can affect on the 
behaviors of the field induced magnetization. First, the modulation was suggested 
as the Ferromagnetic-Ferromagnetic- Antiferromagnetic (F-F-AF) trimerized Heisenberg 
spin-1/2 chains and found that the magnetization curve has a plateau at 1/3 of 
the saturation magnetization value [HI [27] . The F-F-AF trimerized chain has been 
observed experimentally in the compound 3CuCl2-'2dx [28]. During the recent years, 
the trimerized Heisenberg chains have been studied in details [291 EOl [31] . Using the 
density matrix renormalization group (DMRG) method, the observation of the plateaus 
for chains with different spin (5 = 1/2, 1, 3/2, 2) has been reported [29]. A magnetization 
plateau at m = 1/3 of the saturation value exists at low temperature for both, F-F- 
AF and AF-AF-F trimerized spin-1/2 chains [30]. The spin structure factors are also 
calculated for the trimerized cases and found that the static structure factor dose not 
vary with the external magnetic field at the plateau state |31j . 

The other example of the spin chain with mixed (F-AF) and spatially modulated 
exchange is the compound Cu{3 — Clpy)2{N^)2 where (3 — Clpy) = 3 — Chloroyridine, 
which is known as a typical example of a spin-^ F-F-AF-AF tetrameric spin chain [32] . 
In this system, there is a gap from the singlet ground state to the triplet excited states 
in the absence of a magnetic field. The thermodjTiamic properties of the tetrameric 
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spin chains with alternating F-F-AF-AF exchange interactions have been investigated 
numerically |33l [M]. The existence of a plateau at m = 1/2 of the saturation value 
has been observed. The temperature dependence of the magnetization, susceptibility 
and specific heat are studied to characterize the corresponding phases. In recent paper, 
one of the authors considered a different tetrameric spin chain as AF1-F-AF2-F [35] • By 
means of numerical Lanczos method, nonlinear sigma model and bosonization approach, 
it has been found that a magnetization plateau appears at ^ of the saturation value. 
The effects of the space modulation are reflected in the emergence of a plateau in other 
physical functions such as the F-dimer and the bond dimer order parameters, and the 
pair- wise entanglement. 

In this paper, we consider the spin-i Heisenberg chains with alternating exchange 
supplemented by the additional modulation of the one subset of bonds with period 
three, what gives the hexamer modulation of the spin exchange in the chain. Beside 
the very rich quantum magnetic phase diagram, this model has the ability to exhibit 
two new plateaus theoretically, at values 1/3 and 2/3 of the saturation. We restrict 
our consideration by the case where the exchange on the modulated subset of bonds is 
antiferromagnetic and substantially larger then the exchange on undistorted bonds. In 
the presence of strong magnetic field commensurate with dominant AF exchange, the 
model is mapped onto an effective spin- 1/2 XXZ Heisenberg chain in the presence of both 
longitudinal uniform and spatially-modulated fields. This mapping allows us to use the 
continuum-limit bosonization analysis and study the ground state phase diagram of the 
effective spin-chain model. We show that the very presence of additional modulation 
leads to the dynamical generation of two new energy scale in the system and to the 
appearance of four additional quantum phase transitions in the magnetic ground state 
phase diagram. These transitions manifest themselves most clearly in the presence of 
two new magnetization plateaus at magnetization equal to | and | of the saturation 
value. Also, the width of the new plateaus are estimated using the scaling properties of 
the effective theory. 

We confirm the results obtained within the continuum-limit studies of the effective 
model by using the exact-diagonalization of finite chains and perform an accurate 
simulation at zero temperature using the numerical Lanczos method. There are 
two additional gapped phases and corresponding two magnetization plateaus in the 
magnetization curve of our model at values of magnetization equal to | and | of the 
saturation value. Using the numerical technique, we also study the magnetic field effects 
on the bond-dimer order parameter and string correlation function. Finally, treating the 
non-modulated weak bonds as perturbation, we obtain perfect analytical expressions for 
parameters characterizing width of the two additional plateaus. 

The paper is organized as follows. In the forthcoming section, we introduce our 
alternating space-modulated spin-1/2 model and in the strong modulated coupling limit 
derive an effective spin chain Hamiltonian in section HI. In section IV, we present our 
analytical bosonization results. In section V, the results of a numerical simulation are 
presented. Finally, we discuss and summarize our results in section VI. In the Appendix, 
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Figure 1. (Color online.) Schematic representation of spin chains with the hexameric 
modulation of spin exchange which is considered in the paper, (a) corresponds to the 
hexamer distortion of the exchange pattern "A" type and (b) shows the "B" type. 



the details of calculations considering the used perturbation approach are presented. 

2. The model 

The Hamiltonian of the model under consideration on a periodic chain of N site is 
defined as 

N/2 N/2 N/2 

hJ2{s;,-i + sQ, (1) 

n=l n=l n=l 

where Sn is the spin-1/2 operator on the n-th site, J > and J < denote the 
antiferromagnetic and ferromagnetic couphngs respectively. H is the uniform magnetic 
field and JAF{n) is spatially modulated antiferromagnetic exchange. We restrict 
our consideration by following two types of antiferromagnetic exchange modulation 
corresponding to the hexameric distortion of the exchange pattern: the "A" type 

271 

JAF{n) = J%{1 + 5cos{—n)), (2) 

which corresponds to the hexamer modulation of the spin exchanges given in Fig. [T]^a) 
and the "fi" type 

JafH = J%{1 - -j= sin(yn)), (3) 

which corresponds to the following pattern hexamer modulation of the spin exchanges 
given in Fig. [D^b). The uniform spin-exchange coupling between spins separated by odd 
bonds is considered as | J| <^ J^^. For 6 = and J < the Hamiltonian ([1]) reduces to 
the Hamiltonian of alternating F-AF spin-1/2 chains and for J > to the Hamiltonian 
of alternating AF-AF spin chains in a longitudinal uniform magnetic field. 

3. The effective Hamiltonian 

In the considered limiting case of strong antiferromagnetic exchange on odd bonds 
•Jaf ^ l^^l for magnetic field H ~ Jj^p, one can use the standard procedure 
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to map the model ([T]) onto the effective XXZ modeL The easiest way to obtain the 
effective model, is to start from the limit J = 0, where at = the system reduces to 
the set of disconnected pairs of spins in singlet state, located on odd bonds and coupled 
with strong AF exchange JApin). At H ^ spins on each bond either remain in a 
singlet state, \s), with energy Es{n) = —0.75 JApin) or in one of the triplet states, 
and \t-) with energies Efo{n) = 0.25 Jap{j^) and Et±{n) = 0.25 JapI^) =F H, 
respectively. As the field H increases, the energy of the triplet state \t~^), decreases and 
at i/ ~ JAp{n) forms, together with the singlet state, a doublet of almost degenerate 
low energy state, split from the remaining high energy two triplet states. This allows to 
introduce the effective spin operator r which act on these states as [361 [T6] 



(4) 



The relation between the real spin operator S„ and the pseudo-spin operator t„ in this 
restricted subspace can be easily derived by inspection, 

Q± _ _ C± _ 

'-'2n-l ~ ^2n ~ ' 

Using ([5]), to the first order and up to a constant, we easily obtain the effective 
Hamiltonian 
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where 

htff = H-J%- J/4 , (7) 

and 

hlff{n) = - 6J%cos{27m/3) = -hf cos{2TTn/3), (8) 
in the case of "A" type of modulation and 

hlffin) = sin(-n) = hf sin(2vm/3), (9) 

in the case of "B" type of modulation. 

Thus, the effective Hamiltonian is nothing but the anisotropic XXZ Heisenberg 
chain in an uniform and spatially trimer modulated magnetic fields. The anisotropy 
is A = 1/2 (A = —1/2) in the case of chain with antiferromagnetic (ferromagnetic) 
weak bonds J > (J < 0). It is worth to notice that a similar problem has been 
studied intensively in past years [381 iSHl SO]. At 5 = 0, the effective Hamiltonian 
reduces to the XXZ Heisenberg chain in an uniform longitudinal magnetic field. The 
magnetization curve of this model has only saturation plateau corresponding to the fully 
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polarized chain. At if = and J\p ^ \J\ spins coupled by strong bonds form singlet 
pairs and the singlet ground state of the initial spin-chain system is well described by 
superposition of singlets located on strong bonds with magnetization per bond M = 0. 
In terms of effective r-spin model, this state corresponds to the ferromagnetic order with 
magnetization per-site equal to its negative saturation value m = —1/2. In the opposite 
limit of very strong magnetic field H ^ Jaf^ fully polarized state of the initial chain 
with magnetization per strong bond, M = 1, is represented in terms of effective r-spin 
chain as fully polarized state with magnetization per site m = 1/2. This gives following 
relation between the magnetization per strong bond of the initial spin S* = 1/2 chain, 
(M), and the magnetization per site of the effective r-chain (m): M = m + 1/2. 

At 5 ^ 0, the effective model corresponds to the Heisenberg chain in spatially 
modulated fields with period three. In this case, in agreement with standard theoretical 
predictions, additional magnetization plateaus would appear at values of magnetization 
I and I of the saturation value. Below, in this paper we use the analytical and numerical 
tools to analyze the magnetic phase diagram of the model under consideration and 
determine critical values of magnetic fields corresponding to transitions between sectors 
of the ground state phase diagram characterized by different magnetic behavior. 



4. Analytical results 

4.I. The magnetization onset critical field H~ and the saturation field 

The performed mapping allows to determine critical fields H~ corresponding to the 
onset of magnetization in the system and the saturation field H'^ [3Sj- The easiest way 
to get H~ and H'^ is to perform the Jordan- Wigner transformation [H] which maps 
the problem onto a system of interacting spinless fermions with trimerized modulated 
on-site potential: 

IJI ^'"^ \J\ ^'"^ 

Hsf= ± ^ 5Z {^n^n+l + h.C.^ +— pnpn+1 - {po + Pl{n)) pn , (10) 

n n 

where pn = at.'^ni /^o = ^e// — J/'^ = H — J\p — J/2, pi{n) = h\f^{n) and the sign + 
(— ) corresponds to the J > (J < 0). 

The lowest critical field H~ corresponds to that value of the chemical potential p^ 
for which the band of spinless fermions starts to fill up. In this limit, we can neglect 
the interaction term in Eq. (ITOl) and obtain the model of free particles with three band 
spectrum. Below, in this subsection we consider only the case of "A" type of exchange 
modulation with pi{n) = hljrj{n) given by Eq. ([8]). Generalization of this results for 
case "B" is straightforward. 

In this case three bands of the single particle spectrum are given by 

Ei{k) = - H + J% + JV1 + A2 cos ip{k) , (11) 
E^ik) = - H + J% + JV1 + cos {^{k) + 27r/3) , (12) 
Es{k)= - H + + JV1 + A2 cos {^{k) + 47r/3) , (13) 
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where A = 5J\p/ J, 



,,,1 / cos(3A;) + A3 , , ^, 

^{k) = - arccos / | , (14) 

\ /(i + A2)3 



and — 7r/3 < k < n /3. This gives 



H- = J% + JVI + A2 cos {(fiir/S) + 47r/3) at J > , (15) 
H~ = J% - JVl + A2 cos ^(0) at J < . (16) 

To estimate the critical magnetic field H^, which marks the transition into the 
phase with saturated magnetization, it is useful to make a site-dependent particle-hole 
transformation on the Hamiltonian of Eq. (fTO|) : aj^ — t- d^. Up to constant the new 
Hamiltonian reads 

N/2 J N/2 

Hhole = Th^Y. + ^-C.) + 4 5Z PnpUl - (4 + l^lin)) pi , (17) 



n 



where the hole chemical potential /Xq = — /io + <//2. In terms of holes, H'^ corresponds 
to the chemical potential where band starts to fill up, and one can neglect again the 
interaction term. However, the effect of interaction is now included in the shifted value 
of the chemical potential for holes. After simple transformations, we obtain 

He = JIf + i^ + Jx/l + T^ cos(^(0) at J > , (18) 

J 



H+ = J% + -- JVl + A2 cos ((^(vr/S) + Air/S) at J < . (19) 

The spectrum of the system in the free fermion limit f fTT]) - f|T3|) allows to determine 
two other important parameters which characterize the values of magnetization in the 
magnetization curve of the system which the additional plateaus appear and the values 
of magnetic field which correspond to the center of each plateau. Below we consider 
only the case J > 0, however extension to the case J < is straightforward. 

At 1/3-rd band-filling, all states in the lower band E^{k) are filled and separated 
from the empty at £'2 (A;) band by the energy 2mo = E2{k = — -£3(0). Therefore, the 
first magnetization plateau will appear at magnetization equal to 1/3 of the saturation 
value. The magnetic field at the center of plateau is given by 

Hli = Jaf + ^3(0) + mo. (20) 

Analogically at 2/3-rd band-filling, all states in the lower bands E^{k) and E2{k) are 
filled and separated from the empty at Ei{k) band by the energy 2mo = -E'i(7r/3) — 
£'2(7r/3). Therefore, the second magnetization plateau appears at magnetization equal 
to 2/3 of the saturation value and the magnetic field at the center of plateau is given by 

H^2 = AF + E2{7r/3) + mo. (21) 

Since at finite band-filling the effect interaction between spinless fermions cannot be 
ignored the width of plateau differs from its bare value 2mQ. In the forthcoming section 
we use the continuum limit bosonization treatment of the effective spin-chain model (fTOll 
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to determine parameters characterizing the appearance and scale of the magnetization 
plateaus. 



4-2. Magnetization plateaus: Hf and 



C2 



To determine parameters characterizing the appearance and scale of the magnetization 
plateaus, we use the continuum-limit bosonization treatment of the model ([6]). Following 
the usual procedure in the low energy limit, we bosonize the spin degrees of freedom at 
fixed magnetization m and the interaction term becomes 

A, 



m + 



\j^dx(t){x) H ^ sin (^\/AttK (f){x) + (2m + Vjirn 



-i^/iT/Ke{x) 



1 



sm 



\/ATxK(t){x) + (2m + l)7m 



(22) 
(23) 



where Ai and Bi are non-universal real constants of the order of unity [13] and m is 
the magnetization (per site) of the chain. Here (l){x) and 0{x) are dual bosonic fields, 
dtip = VgdxO, and satisfy the following commutation relations 

[(P{x),9{y)]=zQ{y-x), 

[0(x),^(x)] =V2, (24) 

and K{A,m) is the spin-stiffness parameter for chain with anisotropy A and 
magnetization m. At zero magnetization 

^^'^'°'= 2(l-aLosA) - P^' 

Therefore, at m = 0, -R' = 0.75 for J > and K = 1.5 for J < 0. At the transition 
line into the ferromagnetic phase, where the magnetization reaches its saturation value 
f^sat = 0.5, the spin stiffness parameter takes the universal value K{A,0.5) = 1 [H]. 
Respectively for finite magnetization, at < m < nisat and J > 0, the function K{A, m) 
monotonically increases with enhancing m and reaches its maximum value at saturation 
magnetization, K{A, mgat) = 1, while for J < 0, monotonically decreases with increasing 
m and reaches its minimum value at saturation magnetization K{A, 0.5) = 1. 

Using (122!) and (123|) . in the case of "A" type of exchange modulation, we get the 
following bosonized Hamiltonian 



H 



Bos 



+ 



+ 



2naQ 



sin I 27r(m H — jn 
6 



V AirKcf) ] + sin 



27r(m + -)n + VAtiKc, 
6 



(26) 



while in the case of " B" type of modulation - following 



^Bos 



dx\^[{dxcl^f + {dxef 



B 



+ 



J±_ 

27rao 



cos I 27r(m 



]n 



COS I 27r(m 



]n 



(27) 
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From the bosonized Hamiltonian, one easily gets the commensurate values of 
magnetization m = ±1/6 where a magnetization plateau appears. Away from 
this commensurate values of magnetization, arguments of cosine terms are strongly 
oscillating and in the continuum limit, their contribution is irrelevant. Therefore in 
this case, the model is described by the effective gaussian model, indicating on gapless 
character of excitations and continuously increasing magnetization of the chain with 
increasing magnetic field. 

At commensurate values of magnetization m = ±1/6, the cosine term in ( l26l) and 
(!27|) is not oscillating and therefore the modulated magnetic field term comes into play. 
Up to an irrelevant shift on constant of the bosonic field, the generalized Hamiltonian 
which describes system at commensurate magnetization is given by 

where 

m^ = 5J\p and = 2m^/V3. (29) 

The Hamiltonian ( 128|1 is the standard Hamiltonian for the commensurate- 
incommensurate transition which has been intensively studied in the past using 
bosonization [15] and the Bethe ansatz [16]. Below, we use these results to describe 
the magnetization plateau and the transitions from a gapped (plateau) to gapless 
paramagnetic phases. 

Let us first consider /igj^ = 0. In this case, the continuum theory of the initial 
hexamer spin-chain model in the magnetic field H = Jj^p ± J/4 is given by the quantum 
sine-Gordon (SG) model with a massive term ~ h\cos{v47rK(j)). From the exact 
solution of the SG model [IT], it is known that the excitation spectrum is gapless for 
K >2 and has a gap in the interval < Ji < 2. The exact relation between the soliton 
mass M and the bare mass mo is given by [H] 

M = JC{K) (mo/J)'/^'""^^ , (30) 



where 



C{K)- 



-2K' 



r(i - K/A) 



(31) 



v^r(^) L 2T{K/A) 

Here F is the Gamma function and the spin stiffness parameter K depends on the sign 
of J and magnetization m. 

In the gapped phase, the ground state properties of the system are determined by 
the dominant potential energy term ~ cos(v47r^0) and therefore in the gapped phases, 
the field is pinned in one of the vacua 

(0|V47rA>|0) = (2n± l)7r, (32) 

to ensure the minimum of the energy. 

At h^jrj 7^ (i.e. H ^ J\p ± J/4), the presence of the gradient term in the 
Hamiltonian fl28|) makes it necessary to consider the ground state of the sine-Gordon 
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model in sectors with nonzero topological charge. The effective chemical potential 
~ heff\J^dx(f) tends to change the number of particles in the ground state i.e. to create 
finite and uniform density of solitons. It is clear that the gradient term in fl25]) can be 

eliminated by a gauge transformation — > + h^jj-^J^ x, however this immediately 
implies that the vacuum distribution of the field will be shifted with respect of the 
minima f l32|) . This competition between contributions of the smooth and modulated 
components of the magnetic field is resolved as a continuous phase transition from a 
gapped state at |/ie//l < M to a. gapless (paramagnetic) phase at |^e//| > M, where M 
is the soliton mass [15] . This condition gives two critical values of the magnetic field for 
each plateau H^j = H^j ± Mj. 

In the case of chain with AF exchange on weak bonds ( J > 0) this gives 

= Af ± JCiK,) icaf^'-''^^ , (33) 

= Af + J± JC{K,) {caf^^-""'^ ,z = A,B (34) 

where 

= i\:(l/2, -1/6) , and K2 = K{l/2, +1/6) . (35) 

and ca = I, cb = 2/v^. 

Respectively, in the case of chain with ferromagnetic exchange on weak bonds 
(J < 0) we obtain 

Hti = Af + J± \J\C{K[) (q7)'/^'-^^^ , (36) 
= Af ± \-mK) (q7)'/('-^^^ ,^ = A,B (37) 

where 

K[ = K{-l/2,-l/Q) , and K'^ = K{-l/2, +1/Q) . (38) 
Correspondingly the width of each magnetization plateaus is given by 

Ht, - ^ 2M, ^ ( ^^^(^^') ^"^"^Z^L, .-1.2. (39) 

" \ 2\J\C(K;,) (mj/l J|)" at J < 

To estimate the numerical value of the spin stiffness parameter K at magnetization 
m and anisotropy A we use the following ansatz 

K{A,m)=K{A,0) + ^{l-K{A,0)) . (40) 

rusat 

Here, we assume that with increasing magnetization, the spin stiffness parameter 
monotonically reaches its extremum K = 1 at m = rrisat- This ansatz gives that 
Ki = K2 — 0.87 for J > and Ki = K2 — 1.335 for J < 0. It is straightforward 
to get, that in the antiferromagnetic case the width of the magnetization plateau scales 
as 6^'^^ while in the case of chain with ferromagnetic weak bonds as S^'^^. 

In order to investigate the detailed behavior of the ground state magnetic 
phase diagram and to test the validity of the picture obtained from continuum-limit 
bosonization treatment, below in this paper we present results of numerical calculations 
using the Lanczos method of exact diagonalization for finite chains. 
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Figure 2. (Color online.) Difference between the energy of the first lowest level and 
the ground state (Gap) as a function of the magnetic field, H for chains with exchanges 
(a) jo^p, = S = ^ and J = -1, (b) J% = f, S = ^ and J = 1, (c) J% = 6, 
(5 = i and J -1 , (d) = 6, 5 = i and J = 1 and lengths N = 12, 18 and 24. 



5. Numerical simulation 

A very famous and accurate method in field of numerical simulation is known as the 
Lanczos method. In order to explore the nature of the spectrum and the phase transition, 
we use the Lanczos method to diagonalize numerically the model (II]) with periodic 
boundary conditions. 

In this section, to find the effect of a magnetic field on the ground state phase 
diagram of the model, we calculate the spin gap, the magnetization, the string 
correlation function and the bond-dimer order parameter for finite chains with lengths 
= 12, 18, 24 and periodic boundary conditions. 

Since in a critical field, the energy gap should be closed, the best way to find the 
critical fields is the investigation of the energy gap, which is recognized as the difference 
between the energies of the first exited state and the ground state in finite chains. 
In Fig. ([2]), we have presented results of these calculations on the energy gap for the 
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exchanges parameter corresponding to the "A" and "B" types as 

19 2 
(a) J% = — , 5 = — , J = -1, 



3 19 



1 9 9 
(c) J% = 6, 6 = ^, J=-l, 

id) J% = 6, 5 = ^, J = +l. (41) 

and chain lengths = 12, 18, 24. At if = 0, it is clearly seen that the spectrum of 
the model is gapped. As soon as the magnetic field applies, the energy gap decreases 
linearly with H and vanishes at . This is the first level crossing between the ground 
state energy and the first excited state energy. By further increasing the magnetic field, 
three gapless regions 

H-<H<H-, 

HX<H< H+, (42) 
and three gapped regions 

H-<H< H+, 
H-<H< HX, 

H >Ht. (43) 

are seen in Fig. O Oscillations of the energy gap in gapless regions are the result of level 
crossings in finite size systems. One should note that the energy gap in the last gapped 
region, H > H^, growths linearly with the magnetic field H which is an indication 
of the saturated ferromagnetic phase. Therefore in respect to the non modulated case 
6 = [16], two new gapped regions are created by adding 6. To find the critical fields, 
we have used the phenomenological renormalization group (PRG) method [16]. As an 
example, critical fields for the exchanges = ^,5 = ^,7=— 1 are given as follows: 

= 5.67 ±0.01, i7+ = 6.65 ± 0.01, 

H~ = 5.75 ± 0.01, = 5.91 ± 0.01, 

if," = 6.00 ±0.01, ii+ = 6.62 ± 0.01 (44) 

The first insight into the magnetic order of the ground state of the system 
determined by studying the magnetization process. The magnetization along the field 
axis is defined as 

M^ = -J2(^'\^n\Gs), (45) 

n=l 

where the notation {Gs\...\Gs) represents the ground state expectation value. In Fig. El 
we have plotted the magnetization along the applied magnetic field, M^, vs H for chain 
lengths A^ = 12,18,24 and different exchange parameters corresponding to the "A" 
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Figure 3. (Color online.) Magnetization as a function of the magnetic field, H for 
chains with exchanges (a) Jj^p = S — and J = —1, (b) J^^ ^ T' ^ 
and J = 1, (c) Jj^p = 6, (5 = i and J = — 1 , (d) J^p = 6, (5 = ^ and J = 1 and 
lengths N = 12, 18 and 24. In the insets , the magnetization on site as a function of 
the site number " n" is plotted for a value of the magnetic field in the region of plateau 
of l/3A/,at and 2/3Msat for length N ^ 24. 



and "B" types. As it is clearly seen, in absence of the magnetic field H = 0, the 
magnetization is zero. By increasing the magnetic field, the magnetization remains 
zero up to the first critical field H = H~ . This part of the magnetization is known 
as the zero-plateau. This behavior is in agreement with expectation based on general 
statement that in the gapped phases, the magnetization along the applied field appears 
only at a finite critical value of the magnetic field equal to the spin gap. Besides the 
standard zero and saturation plateaus at H < H~ and H > respectively, two 
additional plateaus are seen at = iMsat and = fM^at, where both of them 
satisfy the Oshikawa-Yamanaka-Affleck (OYA) condition. To check that the middle 
plateaus are not finite size effect, we performed the size scaling |19] of their width and 
found that the size of the plateaus interpolates to finite value when — y oo. As 
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Figure 4. (Color online.) String correlation function as a function of the magnetic 
field, H for chains with exchanges (a) J^p =^ T' ~ ^^^'-^ J = 0^) Jaf = T' 
,5 = ^ and ,/ = 1, (c) jOj, = 6, (5 = i and J = -1 , (d) = 6, 5 = i and J = 1 
and lengths N = 12, 18 and 24. 



it has clearly seen in Fig. |2] and Fig. |3l width of plateaus and the mentioned gapped 
regions grow by increasing the modulation S. In the insets of Fig. [3l the magnetization 
on site, = {Gs\S^\Gs) , as a function of the site number n is plotted for some values 
of the magnetic field corresponding to the plateaus at ^Mgat and ^Mgat and length 
= 24. As we observe, the system shows a well pronounced modulation of the on site 
magnetization, where magnetization on weak modulated-bonds is larger than on strong 
modulated-bonds. This distribution remains almost unchanged within the plateau at 
^Msat for H~ < H < H~ and the plateau at ^Mgat for fields H~ < H < 

By analyzing the results on the energy gap, we found that the spectrum is gapped in 
absence of the magnetic field which is one of the properties of the Haldane phase, since, 
we have considered the very strong antiferromagnetic (AF) exchange on odd bonds, 
the Haldane phase cannot be found and it is expected to be the gapped dimer phase. 
This phase can be recognized from studying the string correlation function. The string 
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correlation function in a chain of length A^, defined only for odd /, is [5U] 

2n+Z+l 

OsUl,N) = -{exp{z7i S^)- (46) 

2n+l 

In particular, we calculated the string correlation function for different finite chain 
lengths = 12, 18, 24. Since the model has the U (1) symmetry in presence of a magnetic 
field, we consider the transverse component of the string correlation function. Fig. H] 
presents the Lanczos results on the string correlation function for the chain with lengths 
= 12,18,24 and different exchange parameters corresponding to the "A" and "B" 
types. As can be seen from this figure, aX H < H~, the string correlation function 
Ostri^i N) is saturated and the hexamer chain system is in the dimer phase. As soon as 
the magnetic field increases from the first critical filed, as expected the string correlation 
function decreases very rapidly and reaches zero in the thermodynamic limit, which 
shows that only the dimer phase remains stable in the presence of a magnetic field less 
than The nonzero values of the string correlation function in the region H > , 
are finite size effects and is completely clear that in the thermodynamic limit A^ — > oo 
will be zero. 

An additional insight into the nature of different phases can be obtained by studying 
the correlation functions. Since there are three kind of the space-modulated exchanges 
in our hexameric chain model, we define the following bond-dimer order parameters; 



(7V/6)-l 

D 



d^^^ — — < Gs\SQn+lSQn+2\Gs > 



N 

n=0 
(7V/6)-l 



df^'^ — — < GslSen+s-^Gn+ilGs >, 

n=0 
g W6)-l 

df^^ = — < GslSgn+s-Sen+elG's >, (47) 

n=0 

where summations are taken over the space-modulated antiferromagnetic bonds. In 
Fig. El we have plotted d^^\ d^^^ and d^^'-' versus magnetic field H for chain of length 
A^ = 24 with different exchange parameters corresponding to the "A" and "B" types. 
As it is seen from this figure, at if < H~ spins on all antiferromagnetic space-modulated 
bonds are in a sing let state d^^' = d^^^ = d^^^ ~ -0.75, while at H > H+, the bond- 
dimer order parameter dr, is equal to the saturation value d^^^ = d^^'^ = d^^^ ~ 0.25 
and the ferromagnetic long-range order along the magnetic field axis is present. However, 
in the considered case of strong antiferromagnetic exchanges [Jaf ^ l^^l) ^^"^ high 
critical fields, quantum fluctuations are substantially suppressed and calculated averages 
of spin correlations are very close to their nominal values. 

For intermediate values of the magnetic field, at H~ < H < the data presented 
in Fig. |5] gives us a possibility to trace the mechanism of singlet-pair melting with 
increasing the magnetic field. As it follows from Fig. |5l at values of the magnetic field 
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Figure 5. (Color online.) The bond-dimer order parameter as a function of the 
magnetic field, H for chains with exchanges (a) Jj^p — 'Y^ ^ ^ T§ ^^'^ ^ 
(b) jo^ = ^, 5 = A and J = 1, (c) jo^ = 6, (5 = i and J -1 , (d) jo^ = 6, (5 = i 
and J = 1 and length N = 24. 



slightly above H~ spin singlet pairs start to melt in all modulated bonds simultaneously 
and almost with the same intensity. By further increasing of H and for fields H > H~ , 
melting of weak modulated bonds gets more intensive, however aX H = H~ the process 
of melting stops and the bond-dimer order parameter remains constant up to the critical 
field H = As it is seen in Fig. [5t^a) and (b) for the "A" type of the modulation, 
in the 1/3-plateau state, weak and strong modulated bonds manifest on-site singlet 
features with dimerization values ~ —0.25 and ~ —0.70, respectively. In contrast, for 
the "B" type of the modulation (Fig. [5](c) and (d)), weak modulated bonds are almost 
polarized with dimerization value ^ 0.18, while intermediate and strong bonds manifest 
strong on-site singlet features with value ~ —0.70. By more increasing the magnetic 
field and for fields H > iJ^, the melting of the weak bonds happens intensively and 
aX H = H~ the process of melting stops and remains stable up to the critical field 
if^. It is clearly seen, in the 2/3-plateau state, weak bonds for "A" type and weak and 
intermediate modulated bonds for "B" type are polarized (^ 0.22), but strong bonds 
still manifest strong on-site singlet features (^ —0.65). Finally, for if+ strong bonds 
start to melt more intensively while the polarization of weak bonds increases slowly and 
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at H = all weak and strong bonds achieve an identical, almost fully polarized state. 

6. Conclusion 

In this paper, we have studied the effect of additional modulation of strong 
antiferromagnetic bonds on the ground state magnetic phase diagram of the ID spin-1/2 
chain with alternating exchange. In particular, we focus our studies on the chain with 
hexamer modulation, where the spin exchange on even bonds is the same, while the 
strong antiferromagnetic exchange on odd bonds is modulated with period three. 

In the limit where the odd couplings are dominant, we mapped the model to an 
effective XXZ Heisenberg chain with anisotropy A in an effective uniform and spatially 
trimer modulated magnetic field. The anisotropy parameter is A = ^ in the case of 
antiferromagnetic exchange on even bonds and A = — | in the case of ferromagnetic 
exchange on even bonds. Using the continuum-limit bosonization treatment of the 
effective spin-chain model, we have shown that the additional modulation of the strong 
bonds with period three and amplitude ~ 5 leads to generation of two gaps in the 
excitation spectrum of the system at magnetization equal to the 1/3 and 2/3 of its 
saturation value. As a result of this new energy scale formation, the magnetization 
curve of the system M{H) exhibits two plateaus at M = ^Mgat and M = ^Mgat- The 
width of the plateaus, is proportional to the excitation gap and scales as 6'^, where 
critical exponent i> = 1.13 ± 0.01 in the case of a AF exchange on even bonds and 
u = 1.50 ± 0.01 in the case of ferromagnetic exchange on even bonds. 

For complete description of the model, we supplement our analysis by a very 
accurate numerical simulation. Using the Lanczos method of numerical diagonalizations 
for chains up to = 24, we have studied the effects of an external magnetic field on the 
ground state properties of the system. In the first part of the numerical experiment, we 
have focused on the energy gap of the system. Our results showed that, in respect to the 
non-modulated chain, two new gapped regions create by adding the modulation. The 
widths of the mentioned gapped regions grow by increasing the parameter of modulation 
6. In the second part of the numerical experiment, we have studied the magnetization 
process. We provided a clear picture of the magnetization which showed that two 
magnetization plateaus appear at values ^Mgat and ^Mgat in the new gapped regions. 
To find additional insight into the nature of different phases, we also calculated the on- 
site magnetization, the string correlationfunction and the bond-dimer order parameters. 
The on-site magnetization showed a microscopic picture of the direction of spins on 
different sites, when system is in the new gapped regions. On the other hand, by 
studying the string correlation function, we found that in the absent of the magnetic 
field, the suggested alternating chain is in the dimer phase and this phase remains stable 
in the presence of an external magnetic field up to the first critical field. Finally, we 
studied the effect of the magnetic filed on the ground state phase diagram of the model, 
by means of the perturbation approach. Using perturbation theory, we provided the 
analytical results for the critical fields that these results were in well agreement with 
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the obtained numerical experiment results. 
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Appendix 

Perturbation results 

In this section, we study the effect of the magnetic filed on the ground state phase 
diagram of the model, using the perturbation approach. We try to find the analytical 
results for the critical fields. The behavior of the model in the limit of strong 
couplings on the odd bonds JAF{n) ^ J is interested. For this aim, it is convenient to 
rewrite the Hamiltonian ([1]) in the form 

n =no + v 

N/2 N 
An 

n=l n=l 

N/2 

V = J^S2n.S2n+l. (48) 
n=l 

The unperturbed part, T/o, is the Hamiltonian of N/2 non- interacting pairs of spins. 
The eigenstate of the unperturbed Hamiltonian is written as a product of pair states. 
By solving eigenvalue equation of an individual pair, one can easily find the eigenstates 
as mentioned in Sec. HI. Let us start with the case of 5 = 0. Since the ground state 
energy of a distinct pair is two fold degenerate at if = Jaf^'^) , the ground state 
energy of unperturbed Hamiltonian "Ho is 2^/^ times degenerate [36] . The perturbation 
V splits this degeneracy. By applying the first order and second order perturbation 
theory for finite chains with periodic boundary conditions and generalize results to the 
thermodynamic limit, one can find two critical fields. In the case modulated chain, 
(5 7^ 0, there are different kinds of bonds: strong and weak. In this case, by increasing 
the magnetic field, first weak bonds melt and go to the triplet state in respect to the 
strong bonds. Therefore, it is naturally to find four additional critical fields in respect 
to the non-modulated case, 6 = 0. 

The determined critical fields for Hexameric chain of the "A" type with J < and 
by using the perturbation approach are 
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Table 1. The critical fields which are obtained by the perturbation theory and the 
corresponding values obtained by the numerical experiment in the Hexameric chain of 
the "A" type and J < 0, the "A" type and J > 0, the "B" type and J < and the 
"B" type and J > 0. 



Critical fields Perturbation Results Numerical Results 



'A" type and J < 



H- 5.62 5.67 ± 0.01 

H- 5.75 5.75 ± 0.01 

H+ 5.87 5.91 ± 0.01 

H- 6.00 6.00 ± 0.01 

H+ 6.62 6.62 ± 0.01 

H+ 6.62 6.65 ± 0.01 



"A" type and J > 



H- 5.62 5.61 ±0.01 

H~ 5.75 5.75 ± 0.01 

H+ 6.37 6.40 ± 0.01 

H- 6.50 6.59 ± 0.01 

H+ 7.62 7.61 ± 0.01 

H+ 7.62 7.66 ± 0.01 



'B" type and J < 



H- 4.90 4.90 ±0.01 

H- 4.90 4.96 ± 0.01 

H+ 5.75 5.75 ± 0.01 

H~ 5.75 5.78 ± 0.01 

H+ 6.59 6.60 ± 0.01 

H+ 6.59 6.61 ± 0.01 



'B" type and J > 



H- 4.90 4.87 ±0.01 

H- 4.90 4.92 ± 0.01 

H+ 6.25 6.27 ± 0.01 

6.25 6.33 ± 0.01 

H+ 7.59 7.60 ± 0.01 

H+ 7.59 7.61 ± 0.01 
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(49) 



Generalization of this result for "A" type and J > 0, "B" type and J < and "B" type 
and J > is straightforward. In Table 1, the numerical and perturbation analytical 
results of critical fields were compared for chain of the "A" type and J < 0, the "A" 
type and J > 0, the "B" type and J < and the "B" type and J > 0. The accuracy 
of analytical results are to two significant digits. We emphasize that the perturbation 
results are in well agreement with the obtained numerical experiment results. 
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